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A  nonlinear  syat.ee  exhibited  some  relationship  between 
system  stability  and  controller  off -time.  An  increase  in 
the  off-time  of  the  controller  from  8.5  *  to  78.0  X  resulted 
in  the  nonlinear  system  maintaining  stability  over  50  time 
periods. 


A  nonlinear  system  was  controlled  using  feedback  based 
on  the  unstable  system  parameters.  The  Henon  and  Heiles 
model  was  used  to  provide  an  unstable  periodic  orbit.  The 
linearized  equations  of  motion,  using  Floquet  Theory,  were 
developed  into  variational  equations.  A  controller  was 
derived  baaed  on  a  system  characteristic:  the 
controllability  matrix.  The  controller  gain  was  also 
calculated  based  on  the  controllabilty  matrix. 

Together  with  the  nonlinear  system  the  controller  was 
integrated  for  many  time  periods  using  zero  initial 
displacement.  Stability  was  achieved  for  approximately  50 
time  periods  of  the  system  around  the  zero  displacement 
conditions.  Other  areas  were  investigated,  along  with 
several  other  controller  gains,  without  any  stability  noted. 
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SCALAR  MODAL  CONTROL  OF  A  NONLINEAR  SYSTEM 


USING  TIME  VARYING  GAIN 


CHAPTER  1 


INTRODUCTION 

Control  of  systems  that  are  linear,  especially  conatant 
coefficient  systems,  have  been  explored  in  great  detail. 
Thia  vaat  knowledge  foma  a  baaia  upon  which  one  can  build 
and  expand.  The  expanaion  into  time  periodic  linear  ayateaa 
ia  another  logical  step.  Floquet  Theory  providea  the  tools 
and  analytical  development  needed  to  deel  with  this  next 
step<4:55).  The  nonlinear  system  is  an  extension  of  the 
linear  system  and  the  utility  of  a  linear  system  model 
relies  on  the  validity  of  the  assumption  that  the  model  is 
a  close  approximation  of  the  nonlinear  system.  Linearizing 
a  nonlinear  system  is  done  for  many  reasons:  1.  ease  of 
understanding  the  system,  2.  reduced  complexity  in  dealing 
with  the  system,  3.  tools  are  readily  available  for 
investigation,  and  4.  from  previous  experience  it  is  valid 
< i .e.  it  works  ) . 

One  can  take  a  nonlinear  system,  linearize  it  (details 
later),  and  try  to  control  it.  The  question. 
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"Whether  or 


not  the  control  worked?"  la  closely  tied  to  the  stability 
Information  of  the  system.  Does  the  desired  action  take 
place  even  if  the  setup  wasn't  perfect?  Or  in  perturbation 
theory.  "  Does  the  displacement  from  a  reference  state 
disappear  as  time  goes  on?"  These  are  good  questions,  but 
the  overall  test,  the  teat  that  indicates  success  or 
failure,  is  based  on  the  reactions  of  the  nonlinear  system 
under  the  same  circumstances .  Stability  is  the  grading 
criteria  for  the  control  function.  Stability  is  defined  as 
having  a  tendency  to  return  to  a  desired  state.  A  system 
with  neutral  stability  doesn't  exhibit  any  tendencies  to 
return  to  a  desired  state,  while  an  unstable  system  when 
disturbed  gets  further  and  further  from  the  desired  state. 

AREA  OF  CONCERN 

With  the  background  established,  focusing  can  begin. 
Several  individuals  have  worked  in  the  area  of  cotrolling 
nonlinear  systems  with  time  periodic  coefficients  under  the 
direction  of  Dr.  Calico  and  Dr.  Wiesel  at  the  Air  Force 
Institute  of  Technology.  The  specifics  of  this  report  deal 
with  the  work  started  by  Captain  Liby  in  his  Masters  Thesis 
in  1983(7:1).  Captain  Liby  investigated  a  scalar  modal 
controller  and  how  changing  the  gain  affects  the  region  of 
stability.  In  this  report  an  effort  was  made  to  understand 
and  control  a  nonlinear  time  periodic  system  and  make 
observations  concerning  this  system. 


PROBLEM 


A  nonlinear  ays tea  will  be  chosen  as  a  medium  through 
which  observations  can  be  made.  It  will  be  linearized  and  a 
reference  state  established  for  investigation  using 
perturbation  theory.  The  reference  state  must  provide  a 
single  unstable  mode  upon  which  scalar  control  using  time 
varying  gain  can  be  tested.  Once  the  linear  model  is 
stable,  the  nonlinear  system  with  a  controller  added  will  be 
checked  for  stability. 

OBJECTIVE 

The  objective  therefore  is  to  choose  a  gain  which 
accomplishes  two  sub-goals:  1.  ensure  stability  in  the 
linear  system,  and  2.  stabilizes  the  nonlinear  system  even 
after  many  time  periods.  (Note:  the  second  objective 
requires  the  full  nonlinear  system  plus  the  controller  to  be 
numerically  integrated  over  several  periods.) 


CHAPTER  2 


THEORY 

The  continuation  of  research  in  scalar  modal  control 
relating  the  linear  model  with  the  nonlinear  system  is  the 
main  focus  of  this  report.  Captain  Liby's  thesis  worked  on 
a  scalar  model  controller  of  an  unstable  system  (i.e.  the 
system  had  an  unstable  Poincare  exponent  > .  After 
integrating  the  nonlinear  system  with  the  linear  controller 
Captain  Liby's  system  was  diverging/unstable.  Using  the 
linear  system,  the  controller  moved  the  positive  Poincare 
exponent  into  the  stable  region,  however  the  nonlinear 
system  did  not  exhibit  the  same  tendencies.  The  thrust  of 
my  work  will  be  to  determine  a  gain  function  such  that  the 
nonlinear  system  is  also  stable  after  many  time  periods. 

A  linear  model  was  developed  baaed  on  the  nonlinear 
unstable  periodic  orbit  used  in  Captain  Liby's  thesis.  The 
Henon  and  Heiles  development  was  the  basis  of  this  orbit. 
Floquet  theory  was  used  to  change  from  state  to  modal 
variables  and  scalar  feedback  control  was  used  to  control 
the  single  unstable  mode.  Relationships  relating  the 
nonlinear  equations  to  the  gain  used  for  the  model  feedback 
were  explored.  Time  varying  gain  was  used  and  the 
difficulties  in  finding  the  gain  were  examined. 


HENON  and  HEILES  MODEL 


In  order  to  determine  the  gain  and  solve  the 
instability  problem,  a  test  case  must  be  chosen.  The  Henon 
and  Heiles  model  was  used  for  the  very  reason  that  it  was 
first  developed  because  it  provides  enough  complexity  to 
generate  nontrivial  solutions  yet  still  maintain  simplicity 
for  the  ease  of  understanding  relationships(S:7S> .  The 
potential  function  and  total  energy  in  this  model  are: 

U  <x, y ) 3  l/2<x2^y2*2x2y-2/3y3>  <1> 

E  3  U<x,y)  ♦  1/2 (  <dx/dt >2  ♦  (dy/dt>2  >  <2> 

Changing  to  canonical  units  where  Px  and  Py  are  the 
generalized  momenta  of  x  and  y  respectively,  the  Hamiltonian 
is 

H  3  1/2(x2*  y2  ♦  2x2y  -  2/3y3>  «■  1/2<Px2  ♦  Py2>  <3> 
The  equations  of  motion  are 


dx/dt  * 

d 

H 

/ 

dPx 

*  Px 

<4> 

dPx/dt  3 

- 
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H 

/  dx 

=  -x  -2xy 

<5> 

dy/dt  3 
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H 
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V 

•< 

3  Py 

<6> 

dPy/dt  * 

- 

d 

H 

/  dy 

3  -y  -x2  ♦  y2 

(7) 

Since  there  is  no  time  dependence  in  the  Hamiltonian  then 
the  total  energy  <2)  is  a  constant. 


The  energy  may  take  on  any  value.  Each  value  of  the 
energy  equates  to  a  different  surface  on  which  the  state 
variables  vary.  But  once  the  energy  value  of  the  surface  is 
set.  it  remains  a  constant. 

The  state  variables  used  are  the  position  and  their 
momenta  reapectively . 

*  *  [  x  Px  y  Py  ]T  (B> 

The  equations  of  motion  are  nonlinear  and  can  be 
represented  in  functional  notation  as 

dft/dt  *  f<R)  (9) 

The  time  rate  of  change  of  the  state  is  defined  by  the 
function  f.  When  using  a  controller.  a  reference  state  is 
also  required  to  compare  the  current  state  to  the  desired 
performance.  For  analytical  developments,  a  linear  set  of 
variational  equations  were  developed.  The  assumption  of  a 
small  displacement  from  a  reference  state  must  hold  for  the 
linear  system  (equations  of  variation)  to  be  a  valid 
representation  of  the  nonlinear  system < 11 : 24) . 

FLOQUET  THEORY 

An  unstable  periodic  orbit  of  period  T  was  used  as  a 
reference  state  ft0.  A  small  displacement  aft  was  used  to 
expand  the  f  function  in  a  Taylor's  series  about  the 
reference  state. 


d(*Q  ♦  «*)/ dt  «  £<»0> 


a«d£/d»|«**o  ♦  0(2) 


(10) 


where  0(2)  art  tarns  of  ordsr  2  and  abovo  in  eft. 
jf  Using  aquation  (9)  avaluatad  at  fc0»  aquation  (10)  in  matrix 

fora  is: 

doit/dt  ■  A(t)e«  (11) 

whara  A(t)  *  df/dt|$t«ft0  (12) 

Tha  elaaanta  of  A,  Aj.j  ara  tha  partial  derivativaa  of 
tha  ith  equation  of  motion  with  raspact  to  tha  jth  alamant 
of  the  state  vector.  Equation  (11)  is  a  linear  time- 
periodic  differential  equation*  since  A(t)  has  period  T. 
The  state  transition  matrix  relates  tha  state  at  time  equal 
to  zero  to  any  time: 

*(t)  *  <p(t,0)*(0)  (13) 

The  state  transition  matrix  is  not  periodic*  but  at  t  >  T* 
<P (T »0)  is  called  the  monodromy  matrix.  The  matrix  <p(t*0) 
also  satisfies  the  differential  aquation  (11): 

dq>( t ,0)  /dt  ■  A(t)g(t,0)  (14) 

At  time  equals  zero*  examining  aquation  (13)  reveals  that 
<p(0.0)  is  tha  identity  matrix.  Thus  tha  state  transition 
matrix  can  be  determined  by  solving  aquation  (14)  with  tha 
identity  matrix  as  an  initial  condition. 
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The  Main  result  of  Floquat  ia  <p(t,0)  can  ba  factorad 
uaing  two  natriciea,  F  and  J: 

g( t ,0)  -  F<t>aJtF-l<0)  (15) 

The  Matrix  F(t)  ia  pariodic  with  tha  sane  period  T  aa  the 
original  ayataa.  Tha  J  Matrix  ia  a  conatant  Matrix  which 
when  diagonalized  haa  diagonal  antriea  terned  the  Poincare 
exponenta(2:3) .  Coaplax  axponanta  can  be  repraaanted  aa 
two  real  elenenta  in  Jordan  Noraal  Block  fora.  With  t*T  and 
reMeabering  the  periodicity  of  the  F  Matrix,  o<T,0)  can  be 
calculated  below 

<p(T,0)  *  F<0)aJTF-l(0)  (16) 

where  F<0)  ia  tha  Matrix  of  aiganvactora  of  o(T,0)  (2:5). 

The  eigenvaluea  of  the  Monodroay  Matrix  are  related  to  tha 
Poincare  exponanta  (wj.)  by 

>i  =  a<«iT>  (17) 

Now  F(t)  can  ba  found  by  aubatituting  aquation  (15)  into 
equation  (11)  for  the  atate  tranaition  Matrix  ahown  below. 

dg(t,0)/dt  *  d(F(t)eJtF-l (0) )/dt  (18) 

«  dF(t)/dtCeJtF-l(0)]+CF(t>d(aJt>/dtF-l(0)J  (19) 

Since  F"1(0)  ia  a  conatant,  dF~l(0)/dt  ia  zero.  With 
aoae  rearranging  of  taraa,  tha  tiMe  derivative  of  F(t)  can 


ba  found 


A(t)F(t)eJtF"l(0>  »  dF<t)/dteJtF-l<0>^JF<t>e'JtF-l<0)  (20) 

After  postmultiplying  by  F(0)e~^t,  equation  (20)  becomes 

A(t)F(t>  =  dF(t)/dt  ♦  JF ( t >  (21) 

or 

dF(t)/dt  *  A(t)F (t)  -  F(t > J  (22) 

The  inverse  of  F(t)  can  also  be  found  by  using  the  identity 
relationship  FF"1  *  I  which  with  some  development  becomes 

dF~l(t)/dt  =  -F~l(t)A(t)  ♦  JF~l(t)  (23) 

The  background  is  established  to  introduce  a  change  of 
variables  to  simplify  the  linear  model: 

«(t)  «  F(t)N(t)  (24) 

aubatituting  into  equation  (11) 

dF(t)N(t)/dt  *  N(t)dF(t)/dt  ♦  F(t)dN(t)/dt  (25) 

using  the  latter  relationship,  and  the  right  hand  side  of 
equation  (11)  yields 

A(t)F(t)N(t>*  N(t)dF(t)/dt  ♦  F(t)dN(t)/dt  (26) 

which  when  substituting  in  equation  (22)  yields 

«  N(t) CA(t)F (t)  -  F(t)J)  ♦  F(t)dN(t)/dt  (27) 
Rearranging  and  premultiplying  by  F“l(t>  results  in 


dN(t)/dt 


JN  <  t ) 


(28) 


Thus  the  eigenvector  matrix  raducaa  tha  pariodic  ayataa  into 
a  constant  coafficiant  ayataa (2:7) .  Tha  constant 

coefficient  ayataa  can  ba  aaaily  solved  and  tha  diagonal 
entries  in  tha  J  matrix  (Poincare  exponents)  indeed  give  tha 
stability  information  of  tha  system.  A  positive  diagonal 
element  shows  instability,  were  a  negative  value  shows 
stability.  Since  this  study  is  only  concerned  with  tha 
scalar  control  of  one  mode,  tha  case  of  only  one  mode  with 
a  positive  value  will  be  used.  Feedback  control  is  used, 
but  tha  state  variable  ft,  not  the  model  variable  N  are  used. 
The  state  variables  are  readily  available  and  have  physical 
representations  (  i.e.  velocity). 

CONTROLLER 

Equation  (11)  is  modified  to  include  feedback  control: 

dft/dt  ■  A(t)8(t)  ♦  B(t)u(t>  (29) 

where  B(t)  determines  the  distribution  of  the  control  and 
can  be  a  function  of  time  with  the  same  period  T,  u(t)  * 
G(t)ft(t)  the  feedback  control,  and  G(t>  the  gain  function. 
B(t)  was  chosen  to  affect  only  the  momenta  states. 
Physically  this  relates  to  changing  the  velocity,  whereas 
the  position  states  can  not  be  changed  instantaneously. 
Therefore  B(t)  is  in  the  form 
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*  tO  1  0  1  1  T 

substituting  squation  (24)  into  (29)  rssults  in 

d(F(t)N(t))/dt  *  A(t)F (t)N(t)  ♦  B(t)G(t)F(t)N(t) 
dF (t) /dtN(t ) +F (t)dN(t ) /dt  »  CA ( t ) F (t> *B ( t )G( t )F ( t ) ) N (t> 


(30) 


(31) 

(32) 


rearranging  and  substituting  for  dF(t)/dt,  squation  (22) 

dN(t)/dt  -  tJ  ♦  F-l(t)B(t)G(t)F(t))N(t)  (33) 


Let  g(t) *F"1 (t)B(t)  ba  dsfinad  as  tha  controllability  aatrix 
and  k(t)  *  G(t)F(t)  ba  dafinad  as  tha  naw  gain  aatrix 
Only  tha  unstabla  aoda  will  ba  fadback.  Tha  controllad 
systaa  aquations  taka  tha  fora  balow,  whara  Nc  is  tha  closad 
loop  state. 


dNc/dt  » 


Wl 

0 

*91 

0 

0 

«2 

kg2 

0 

Nc 

0 

0 

kgi*wi  0 

0 

0 

*9n 

wn 

(34) 

Only  tha  ith  alaaant  of  tha  closad  loop  systaa  has  a 
different  Poincare  exponent  froa  tha  open  loop  systea(2:9). 
Tha  unstabla  aoda  aquation  (35)  can  ba  solved  by  tha  use  of 
an  integrating  factor,  yielding  (36): 

dNCi/dt  *  t  wi  ♦  k(t)gi(t>  ] NCi  (35) 


NCi (t) «NCi (0>axp Cq  /  t (wi*k ( t >gi (t) )dt) 


(36) 


Both  9i<t)  and  k(t)  can  ba  periodic  and  expressed  as 
Fourier  series.  The  multiplication  can  be  expanded  as  a 
constant  part  and  a  periodic  part.  The  state  transition 
matrix  is  the  exponential  part  of  equation  (36>;  when 
evaluated  at  the  period  T,  the  periodic  part  goes  to  zero 
<exp<0)=l>  and  stablility  information  is  contained  in  the 
eigenvalue.  Hence  a  new  modal  exponent  is  found. 

Wi'  *  wj  ♦  kC9ic  < 37) 


where  wj/  is  the  new  modal  exponent,  and  c  subacipt 
indicates  the  constant  part  of  gj,  and  k. 


Multiplying 

two 

Fourier 

series  can  be 

resolved 

into 

separate  sine 

and 

cosine 

function  and 

higher 

order 

harmonics.  A  constant  part  will  be  contributed  from  many 
products  of  higher  harmonics,  but  csn  readily  be  reduced 
into  sine  and  cosine  components  with  the  trigonometric  sum 
and  difference  formulas. 

This  is  one  method  for  determining  the  gain,  k(t).  The 
merits  of  this  method  will  be  discussed  later.  Another 
method  for  determining  the  gain  was  developed  because  of  two 
reasons.  First,  the  difficulty  in  multiplying  two  Fourier 
series.  Second,  special  care  must  be  taken  so  the  linear 
assumptions  are  not  violated  even  though  the  linear  system 
is  stable.  Capt  Liby,  in  his  thesis,  chose  a  gain  and 
stablized  the  linear  system,  however  the  nonlinear  system 
was  not  stable  in  all  cases  over  many  time  periods.  Capt 
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Liby  used  a  gain  equal  to  a  constant  tines  the  first 
harmonic  sine  term  (7:16).  By  substituting  this  gain  into 
equation  (36),  he  found  a  new  Poincare  exponent  through 
equation  <37) .  However,  even  though  the  integral  of  the 
remaining  periodic  functions  in  equation  (36)  are  evaluated 
over  one  period  and  hence  zero,  during  the  period  the  values 
can  become  very  large.  When  the  modal  variable  is  large  the 
variation  becomes  large  (7:39).  Hence,  the  small 
displacement  assumption  might  be  violated  by  the  transient 
response  of  the  linear  system.  Once  the  assumption  of 
linearity  is  invalidated,  the  linear  system  is  no  longer  a 
valid  model  of  the  nonlinear  system. 

The  gain  selected  was  developed  by  examining  equation 
<36)  for  the  modal  variable  as  a  function  of  time.  The 
k(t)  function  was  a  discrete,  evenly  spaced  inverse  of  the 
gi  function  multiplied  by  the  desired  constant  so  that  the 
sum  of  kgi  and  wi  is  now  a  negative  (stable)  number. 
Explicitly,  this  drives  the  entire  k(t)gj,(t)  «  constant.  If 
this  can  be  done  successfully,  there  would  not  be  any 
periodic  function  to  cause  the  modal  variable  to  become 
large  during  the  period.  The  variation  will  not  be  large, 
and  the  assumptions  of  linearity  will  remain  valid. The  only 
discrepency  is  the  zero  points  of  the  gi  function.  A  maximum 
value  of  the  k  function  is  equivalent  to  limiting  the 
control  forces  Imposed  on  real  systems.  A  zero  value  can 
be  used  where  a  higher  value  is  required;  this  equates  to  a 
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shutoff  of  ths  control  mechanism  when  the  actual  control 
forces  required  approach  infinity. 

The  particular  value  of  the  constant  (k<t)gi<t>> 
depends  on  several  factors.  Certainly  the  amount  or 
tendency  of  instability  already  present  in  the  system  will 
impose  restrictions  on  the  k<t>gi<t>  value.  Obviously,  if 
during  one  time  period  there  is  not  a  tendency  to  return  to 
the  reference  orbit,  the  deviation  would  grow  in  succeeding 
periods.  This  fact  provides  a  starting  point  for  using 
values  of  the  k(t)gi<t),  it  must  be  at  least  as  negative  as 
the  current  wj,  < Poincare  exponent)  is  positive. 

Another  factor  relates  to  the  off-time  of  the  control 
mechanism.  When  a  zero  point  of  gi<t>  is  reached  during  the 
orbit  the  controller  is  turned  off.  Otherwise  an  infinite 
control  force  would  be  required.  The  futility  of  using  a 
finite  control  input  when  an  infinite  response  is  self- 
evident.  Therefore,  a  trade-off  is  developed  between  the 
percentage  of  off-time  and  the  amount  of  stability 
introduced  by  the  controller  (i.e.  how  big  was  the  constant 
k(t)gi<t)). 

A  final  factor  to  be  considered  is  whether  the  gain  is 
realistic.  A  control  system  that  works  while  requiring 
elaborate  and  extremely  complex  system  operations  may  not  be 
an  overall  improvement  but  a  hinderence.  This  certainly  is  a 
lessor  criterion  but  should  be  kept  in  mind  when  developing 
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Once  the  gain  function  ia  calculated  the  next  atep  in 
the  procedure  ia  incorporating  the  gain  in  the  nonlinear 
systea.  Refer  back  to  equation  (9).  It  now  can  be 
rewritten  with  a  control  function  added. 

dR/dt  *  f (St)  ♦  Control  (3B) 

where  the  Control  ia  a  n  x  1  matrix.  The  composition  of  the 
control  function  ia  determined  by  going  through  the 
transformation  of  (1)  feeding  back  the  atate  variables.  (2) 
converting  to  modal  variables,  <3)  applying  scalar  control 
of  the  unstable  mode,  and  finally  (4)  multiplying  by  the 
gain  function.  Thia  can  be  outlined  by  rearranging 
equation  (24)  in  terma  of  displacement  from  a  reference 
orbit. 


N  ( t ) 

=  F-l(t)St(t) 

(24) 

R(t> 

*  tki(t)  0 

0 

0 

J 

(39) 

where  R(t)  assures  the  gain  of  only  the  unstable  mode. 
B ( t ) ,  equation  (30)  is  uaed  to  properly  distribute  the 
controller  in  the  proper  order, 

B  *  C  0  1  0  1  (30) 

combining  all  terma. 


Control  =  B(t)R(t)F“l (t ) C  x(t> 


xQ(t)  ] 


(40) 


which  reduces  to: 


Control  *  kiC£i"l<t>  IC  x(t)-Xo<t>] 

where  fi~l<t>  is  the  row  of  F"l<t>  associated 
unstable  mode,  ki  is  the  gain,  and  x(t>-x0(t) 
displacement  from  the  reference  orbit  at  any  time. 


CHAPTER  3 


RESULTS 


Using  the  energy  function  (2)  of  the  Henon  end  Heiies 
problem,  an  unstable  periodic  orbit  was  calculated.  Initial 
conditions  of  energy*  0.125  ,  Py<0)=0,  x<0)=0  were  used  with 
an  estimate  for  y(0)  and  the  period  T.  A  predictor- 
corrector  numerical  integrator  was  used  to  match  end  points 
on  the  orbit.  The  period  of  the  orbit  and  initial  state 
values  are  found  in  table  1 . 


ORBIT  VALUES 


T  =  32.040441196 
E  *  .124999316 
X(0)  =  .000000134 
PX<0>  =  .499704388 
Y  (0>  *  -.017054328 
PY (0)  =  .000001441 


TABLE  1. 


The  values  for  the  reference  orbit  for  each  element  in 
the  state  x,  Px,  y,  Py  are  shown  in  Fig  1  -  Fig  4 
respectively.  These  figures  show  50  points  evenly  spaced 
over  the  period.  All  other  figures  that  are  used  to  show 


stability  information  are  deviations  from  these  values 
throughout  the  period. 
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Fig  4.  Reference  Orbit  -  state  variable  Py 
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In  conjunction  with  finding  the  periodic  orbit,  the 
state  transition  matrix  (Id)  was  caiculatsd  at  sach 
intsgration  stsp  and  ths  aonodroay  matrix,  with  ths 
associatsd  stability  information  was  availabla  at  t  *  T. 
Ths  sigsnvaluss  (17)  ars  ralatsd  to  ths  Poincars  sxponsnta 
and  rsvsalsd  a  singls  unstabls  mods.  Ths  sigsnvsctor 
associatsd  with  ths  positivs  Poincars  sxponsnt  is  a  column 
of  ths  F(0)  matrix.  Ths  F(0)  matrix  was  invsrtsd  and  ths 
row  aquation  (42)  corrssponding  to  ths  unstabls  sigsnvslus 
was  ussd  as  an  initial  condition  for  aquation  (23) .  Ths  F~ 
1 ( t )  row  was  intsgratsd  for  ons  psriod  T.  Sines  F-1(t>  has 
ths  sama  psriod  as  ths  linsar  systsm,  ones  it  is  caiculatsd 
for  ons  psriod  it  is  known  for  all  tims.  Ths  row  of  F"1(0) 
is  shown  balow  where  fi“l(0)  is  ths  row  of  F_1  associatsd 
with  wi. 

wi  =  0.083668212  (42) 

fi"l(0)=  C.013240546  .824668173  .417510018  -.381684565)  (43) 

One  must  remember  the  controllability  matrix  g(t)  is 
defined  below: 

g(t )  *  F-l(t)B(t)  (44) 

The  controllability  function  matrix  ussd  in  aquation 
(35)  reduces  aquation  (36)  to  control  of  ths  unstable  mods. 
Ths  gi  in  aquation  (36)  can  be  found  by 


9i<t>  * 


(45) 


v 


where  9i <t)  is  a  1  x  1  matrix,  i.a.  scalar  modal  control. 
Hanes, 

9i<t)  *  fi2_1<t>  «■  fi4-l<t>  (46) 

Tha  scalar  gi<t>  is  tha  sum  of  the  second  and  fourth  element 
of  the  row  of  F~l(t)  corresponding  to  the  positive  Poincare 
exponent . 

In  all  of  the  above  computations  periodic  functions 
were  represented  by  Fourier  series  (1:109>.  The  Fourier 
series  for  gi<t),  first  six  terms,  are  shown  in  Table  2. 


gi<t)  Fourier  series 
coa<nw0t)  sin(nw0t> 


a0  • 0000077 


«1 

. 3300262 

bl 

.4913984 

«2 

.1247512 

b2 

-.1374558 

®3 

.0000185 

b3 

-.0000165 

«4 

-.5316386 

b4 

-.1825139 

«5 

. 3743597 

bs 

-.4790797 

«6 

- . 0000052 

b6 

-.0000148 

TABLE  2. 

Once  the  controllability  matrix  is  calculated  the  linear 
system  can  be  stabilized  with  a  judicious  choice  of  gain. 
Recall  equation  <36>,  and  if  the  sum  wi  ♦  k(t)gi(t>  is  less 
than  zero  then  the  controlled  modal  variable  Nci<t> 
displacement  decays  with  time.  This  implies  the  system 
returns  to  the  reference  state;  stability  is  achieved.  For 
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this  specific  csss  wi  (42)  is  slightly  positivs.  Zf  k(t) 
and  9i<t>  are  rvpresentad  by  fouriar  aeries,  tha 
multiplication  of  tha  two  in  ganaral  raquiraa  multiplying 
two  infinita  sarias.  Zf  n  tarns  ara  usad  for  each,  k(t) 
thaoratical ly  could  ba  calculatad  from  2n  ♦  1  linaar 
aquations.  Thara  would  ba  n  cosina  and  sina  aquationa  and 
ona  constant  part  aquation  (37).  Tha  idaal  rastriction  of  a 
zero  coafficiant  for  tha  first  n  tarns  was  axplorad.  This 
would  provide  for  the  first  n  higher  harmonics  to  be  zero. 
Specifically.  28  terms  for  a  total  of  57  aquations  ware  usad 
to  determine  a  k(t).  The  new  exponent  w(t)  was  plotted  and 
showed  non-zero  higher  harmonics.  high  ocillations.  The 
reason  for  the  poor  outcome  is  obvious  when  one  rearranges 
tha  variables  in  exponent  of  aquation  (38) . 

k(t)  *  (  wi'  -  wi  )/  gi(t>  (47) 

where  wi  is  the  old  Poincare  exponent  and  wi '  is  tha  desired 
Poincare  exponent  (negative  for  stability).  When  a  zero  of 
gi(t)  is  reached  k(t)  is  undefined.  A  representation  of  k(t) 
shows  this  descrepency  as  g(t)  approaches  zero  from  plus  to 
minus,  k(t)  changes  from  a  large  negative  value  to  a  large 
positive  value  over  a  short  time  change.  The  graph  of  gi(t) 
is  shown  in  Fig.  5  and  the  k(t)  determined  in  Fig.  6.  Zf 
one  looks  closely  the  spikes  on  the  k(t)  graph  correspond  to 
tha  zeroes  on  tha  gi(t)  graph. 


Controllability  scalar  for  unstable  sods 


An  attempt  was  mada  to  aliainata  all  of  tha  zaroaa  of 
howavar,  this  is  ispossibla  bacausa  of  the  coaposltion  of 


gi<t).  Since  gi<t>  is  a  linear  combination  of  f2"*<t>  and 
f4~l<t)  one  could  dataraina  a  ratio  of  tha  two  to  limit  tha 
zeroes  to  only  two.  But  tha  limit  is  two  bacausa  tha 
functions  sine  and  cosine,  have  two  zeroes  for  each  period. 
With  this  in  mind,  gi<t)  was  examined  and  two  restrictions 
to  finding  kCt)  ware  established:  (1)  a  maximum  value  of  10 
units  was  sat  for  k<t>,  and  (2)  whan  a  higher  value  of  k(t) 
is  determined  from  equation  (47)  then  k(t)  was  set  to  zero. 
This  is  just  a  simple  restriction  so  an  infinite  input  is 
not  required  from  the  controller.  Following  these 
restrictions  and  using  1000  discrete  evenly  spaced  values 
between  zero  end  the  period,  k<t>  is  shown  in  Fig.  7. 


modal 

value 

exponent, 

function . 

from  equation 

(36) ,  is  now 

reduced 

Wi' 

*  Wi 

unstable , 

k  *  zero 

(43) 

Wi' 

=  wi  -  kgi 

stable,  k 

non-zero 

(49) 

The  percentage  of  the  time  out  of  one  period  when  k  is 
zero  is  inversely  proportional  to  the  maximum  value  of  k(t). 
The  desired  value,  wi'  can  compensate  for  this  off -time. 
The  restrictions  impossed  on  the  gain  function  simplify  the 
equation  (36)  for  the  displacement  of  tha  modal  control 


variable. 


It  now  is  either  equation  (50)  or  (51). 

Nci(t)  =  NCi (0)ewi ' t  (50) 

Nci(t)  =  Nci(0)ewit  ,k  =  0  (51) 

From  equation  (50)  and  (51)  one  can  clearly  ahow  that 
the  displacement  is  either  decaying  or  slightly  growing.  A 
tradeoff  exists  between  the  amount  of  decay  and  growth  based 
on  relative  magnitude  of  wj.'  and  wj.,  which  most  of  the  time 
decays.  The  choice  of  gain  may  seem  complicated  but  the 
goal  la  to  stabilize  the  nonlinear  system  as  well  as  the 
linear  solution. 

TESTING 

Eight  different  gain  configurations  were  tested  and  are 
shown  in  Table  3.  The  nonlinear  system  was  numerically 
integrated  for  several  periods,  each  period  consisting  of 
1000  evenly  spaced  time  steps.  Prior  to  the  first  timestep 
the  displacement  of  the  current  state  from  the  reference 
state  is  recorded.  Thereafter  for  each  group  of  20 
timesteps  for  a  total  of  50  deviations  are  recorded  for 
every  period  integrated.  For  example,  if  the  nonlinear 
system  is  integrated  over  18  periods  and  50  deviations  are 
recorded  per  period  then  there  will  be  a  total  of  900  state 
displacements.  These  900  points  can  be  plotted  for  each 
variable  of  the  state  to  provide  a  graphic  representation  of 


the  trend 


Patterns  can  easily  be  seen  and  the  controller 


effectiveness  can  also  be  evaluated. 


TESTING  CASE  DATA 


MAGNITUDE 

k(t)g(t) 

MODAL  EXPONENTS 

CASE 

k ( t) max 

g(t)  CUTOFF 

(constant) 

off 

on 

1 

0 

_ 

0 

♦ .08367 

— 

2 

10 

0.02641 

-  .2641 

♦ . 08367 

-  .18045 

3 

11 

0.03601 

-  .3961 

♦ .08367 

-  .31252 

4 

10 

0 . 05682 

-  .5682 

♦.08367 

-  .44458 

5 

1 

0.10000 

-  .1000 

♦  .08367 

-  .01633 

6 

10 

0.10000 

-1.0000 

♦.08367 

-  .91633 

7 

10 

0.20000 

-2.0000 

♦ .08367 

-1.91633 

8 

2 

1 . 00000 

-2.0000 

♦  .08367 

-1.91633 

TABLE  3. 

Each  case  is  comprised  of  four  main  components: 
k(t)max,  g<t)  cutoff,  k(t)g(t)  constant  ,  and  modal 
exponent.  The  k<t)max  is  the  maximum  magnitude  of  the  gain 
function.  An  arbitrary  cutoff  of  around  10  was  used  to 
simulate  an  acceleration  limitation. 

The  g(t)  cutoff  determines  the  on/off-time.  During  the 
gain  function  calculation  in  equation  (47)  when  the 
magnitude  of  g(t)  falls  below  the  g(t)  cutoff  value,  the 
k(t)  gain  function  is  assigned  the  value  of  zero.  The  lower 
the  g(t)  cutoff  the  less  off-time.  Case  6  with  a  g(t) 
cutoff  of  0.1  has  only  6.5k  off-time.  Whereas  case  8,  g(t> 
cutoff  of  1.0,  has  an  off-time  of  78.0k 

The  constant  value  ,  k(t)g(t),  determines  the  tendency 
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for  decay  of  the  displacement  of  the  modal  variables.  (see 
equations  (48)  and  (SO)  )  The  constant  value  can  be  used  to 
compensate  for  increased  off -time. 

The  modal  exponent,  as  defined  by  equations  (48)  and 
(49),  correspond  to  the  off  and  on  time  of  the  controller. 
The  modal  exponent  values  range  from  a  low  of  -.016  to  a 
high  of  -1.9  in  terms  stability  when  the  controller  is  on. 

Case  1  is  with  zero  gain,  so  the  deviation  of  the 
nonlinear  system  to  zero  initial  displacement  from  the 
reference  state  can  be  examined.  The  numerical  roundoff  and 
approximations  excite  the  unstable  mode.  Fig  8.  shows  the 
deviations  over  27  time  periods  of  the  first  state  variable. 
All  four  state  variables  exhibited  the  same  characteristics 
so  only  the  X  state  variable  deviation  is  shown. 

Cases  2-4  were  used  to  compare  results  obtained  by  Capt 
Liby  using  a  different  controller.  One  can  not  compare 
results  directly  since  Capt  Liby  used  a  sinusoidal 
controller  and  all  above  cases  tested  have  on/off 
characteristics.  Capt  Liby  found  areas  of  stability  and 
these  were  chosen  as  starting  points.  Several  other  points 
in  the  state  space  including  zero  displacement  from  the 
reference  orbit  were  investigated.  Deviations  on  the  order 
of  10~1  were  produced  and  no  stable  areas  could  be  found. 

Cases  5-6  were  chosen  because  Capt  Liby  could  not  find 
any  stable  areas  with  comparable  modal  exponents.  However, 
although  similar  deviations  were  produced  as  above,  no 


r 


o 


© 


t 


02; 


Fig  9.  X  state  deviation  Case  5,  Zero  Displacement 
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stable  areas  could  be  found.  Fig  9.  shows  a  typical  trend 
of  cases  2-6.  Again  only  the  X  state  variable  is  shown. 
Specifically,  Fig  9.  is  case  5  integrated  for  14  tine 
periods  with  no  initial  displacement  from  the  reference 
state. 

Case  7  and  8  are  an  offshoot  of  the  previous  cases. 
The  item  that  differentiates  these  case  from  the  others  is 
the  amount  of  off-time.  The  highest  off-time  in  the 
previous  cases  was  case  6  with  6.5*  off -time.  Case  8  has  an 
off-time  of  78.0*,  or  only  22.0*  on  time.  Both  cases  have 
the  same  k(t)g(t)  constant,  but  differ  in  the  g  cutoff. 

Case  7  is  more  similar  to  the  previous  cases.  Several 
areas  were  investigated  including  zero  initial  displacement 
from  the  reference  state  and  no  areas  of  stability  were 
found . 

Case  8  was  more  promising.  Fig  10.  portrays  the  gain 
function  k(t>  for  case  8.  When  compared  to  Fig  7.  the 
increase  in  off -time  is  dramatic.  Several  areas  were 
investigated  and  stability  was  found.  Fig  11.  shows  zero 
initial  displacement  integrated  for  27  time  periods  with 
decreasing  deviations  on  the  order  10~2.  Again  only  the  X 
state  deviations  are  shown  as  the  other  state  variable 
characteristics  are  similar.  One  can  compare  Fig  11.  to  Fig 
8.  for  evaluating  the  case  8  controller  effectiveness, 
(note  different  scale)  Fig  8.  scale  is  67  times  bigger  than 


In  Fig  11.  the  deviations  do  not  disappear  after 
tapering  off.  To  ensure  the  stability  of  case  8,  the  zero 
initial  conditions  were  integrated  for  77  time  periods, 
shown  in  Fig  12.  After  roughly  50  time  periods  with  the 
deviations  on  the  order  of  10-2  the  deviations  suddenly  grow 
and  the  system  exhibits  unstable  behavior.  This  unstable 
behavior  is  assumed  to  be  a  numerical  mismatching  due  to  a 
possible  synchronization  problem  between  the  actual  period 
of  the  nonlinear  system  and  the  estimated  period  used  in  the 
integration.  With  an  incorrect  period  the  controller  would 
not  function  at  the  proper  time  causing  increased 
displacements  instead  of  reduced  displacements.  However, 
the  overall  significance  of  Fig  12.  is  the  depiciton  of  the 
stable  deviations  over  many  time  periods  (roughly  50> . 


CHAPTER  4 


CONCLUSIONS 

The  title  of  this  report  implied  a  goal  to  be  reached 
at  the  end  of  this  research  endeavor;  control  of  a 
nonlinear  systes.  This  goal  was  to  be  reached  by  solving 
the  stated  problem  and  accomplishing  the  objective  (see  p4> . 
The  measure  of  success  depends  on  success  in  reaching  the 
objective  and  problem  solution.  Therefore,  objective  and 
problem  solution  will  be  exsmined  first. 

The  objective  was  to  choose  a  gain  which  accomplished 
two  sub-goals:  1.  ensure  stability  in  the  linear  system,  and 
2.  stabilize  the  nonlinear  system.  In  a  system  with  a 
single  unstable  mode,  current  theory  allows  linear  systems 
to  be  easily  stabilized.  Once  the  Poincare  exponent  is 
found  a  scalar  constant  gain  can  be  used  and  the  new  system 
is  stable  (refer  to  equation  (37)  pl3) .  However,  with  the 
nonlinear  system  also  dependent  on  choice  of  gain  for 
stability  some  gains  that  stabilize  the  linear  system  are 
rejected  on  the  basis  of  not  accomplishing  the  second  sub¬ 
goal.  Case  2  through  7  fall  into  this  category.  Before 
final  judgment  is  pronounced  on  case  S,  permit  an 
examination  of  the  problem  and  methodology  of  the  solution. 

The  methodology  for  the  problem  solution  relied  on  the 
Henon  and  Heiles  problem  to  act  as  a  test  basis  (5:73).  The 


Henon  and  Hailes  potential  function  was  uaed  for  its  level 
of  complexity  yet  ease  of  use.  The  linearized  equations  of 
motion,  using  Floguet  theory,  were  developed  into  equations 
of  variation.  A  controller  was  derived  based  on  system 
characteristics:  the  controllability  matrix,  g<t).  A  first 
attempt,  baaed  on  multiplying  two  Fourier  series  together  to 
obtain  the  gain  function  k(t>,  was  rejected  due  to  high 
oscillations  causing  further  instability  in  the  nonlinear 
system.  A  second  approach  toward  gain  selection  was  more 
direct  <refer  to  equation  (47)  p24)  Inversion  of  g(t).  and 
multiplying  by  a  constant.  This  is  the  basis  of  cases  2-6. 
A  refinement  of  this  approach  resulted  in  the  last  attempted 
gain  selection,  cases  7  and  8.  As  previously  stated,  each 
gain  technique  was  uaed  in  the  nonlinear  system  via  a 
controller  and  checked  for  stability. 

The  methodology  was  straightforward  and  directly 
attacked  the  problem.  A  nonlinear  system  with  a  single 
unstable  mode  was  controlled  using  a  modal  controller.  A 
technique  for  gain  selection  was  developed 

k<t)  *  <  wi'  -  wj.)/  gj.<t>  )  (47) 

where  k(t)  is  the  gain,  w*  =  Poincare  exponent  /  old  modal 
exponent,  w*'  =  desired  exponent  /  new  modal  exponent  and 
gi<t)  part  of  the  controllability  matrix.  In  this  case  the 
cutoff  value  of  gi<t>  was  also  recognized  as  significant. 
The  cutoff  value  of  gi(t>  determined  the  on/off-time  of  the 


controller.  This  ranged  from  6.5  X  in  case  6  to  78.0  x  in 
case  8.  Case  8  was  the  only  case  that  stabilized  the 
nonlinear  system. 

Was  the  uniqueness  of  case  8  the  ingredient  that  made 
it  the  only  succesful  case?  Clearly  the  off-time  is  a 
marked  increase  over  the  other  cases  (refer  to  Fig  10.  p36> . 
The  results  of  using  this  gain  function,  active  for  only 
22.0  X  of  the  time,  reduces  the  deviations  in  Fig  8.  by  two 
orders  of  magnitude  to  Fig  11.  <p32  &  p37). 

With  this  specific  problem  solved  and  the  objective 
accomplished  now  let  us  go  back  to  measuring  the  success  of 
reaching  the  primary  goal.  Yes,  a  nonlinear  system  was 
controlled.  But  only  in  a  limited  manner.  This  report  does 
not  profess  to  having  a  methodology  to  control  all  other 
nonlinear  sytems,  yet  it  does  make  these  statements: 

1.  It  reaffirms  the  usefulness  of  the  Henon  and  Heiles 
paper  as  a  test  bed  for  experimentation. 

2.  Floquet  theory  is  a  useful  tool  by  which  modal 
controllers  can  be  implemented. 

3.  A  new  gain  selection  technique  based  on  the  system 
controllability  matrix  may  provide  future  utility. 

4.  In  this  specific  case,  a  relationship  between  the 
unstable  Poincare  exponent  and  the  gi<t)  cutoff  value 
or  on/off -time  was  uncovered. 

5.  The  above  relationship  impacted  heavily  on  the 
effectiveness  of  the  controller  and  the  ability  to 
stabilize  the  nonlinear  system. 


From  these  conclusions  several  recommendations  arise. 


RECOMMENDATIONS 


First,  before  any  continued  research  is  conducted  a 
small  investigation  into  whether  the  results  of  this  report 
can  be  generalized  for  the  area  of  modal  control  of  I 

nonlinear  systems  should  be  undertaken.  Specifically,  (1) 
whether  the  reference  state,  with  its  correspondind  unstable 
mode,  can  be  classified  as  a  special  case,  and  (2)  whether 
the  on/off -time  characteristics  of  the  controller  depends  on 
this  classification. 

Next,  if  the  results  were  not  of  a  special  case  and 
could  possibly  be  generalized  then  a  major  decision  must  be 
made.  Either  use  the  Henon  and  Heiles  paper  or  use  an 
unstable  nonlinear  system  with  physical  representations. 

Once  this  decision  is  made,  I  suggest,  a  trade-off  study 
using  the  unstable  Poincare  exponent  and  the  on/off-time 
<gi < t )  cutoff  value)  as  parameters. 

Although  research  is  being  conducted  involving  two  or 
more  unstable  modes  ,  I  believe  investigation  of  the  trade¬ 
off  study  should  not  be  neglected.  It  might  provide  some 
insight  or  other  help  in  conducting  the  multiple  unstable 
mode  research. 

Finally,  in  a  side  note,  I  petition  anyone  in  a 
position  of  authority  with  a  technical  problem  to  solicit 
engineering  graduate  schools  for  research.  The  masters 
thesis  is  an  important  part  of  the  learning  process  in 
higher  education  and  may  provide  some  reciprocal  benefits. 
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